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ABSTRACT 

Using Spitzer IRAC mid-infrared imaging from the Great Observatories Origins Deep Survey, we study z&so- 
dropout sources in the Hubble Ultra Deep Field. After carefully removing contaminating flux from foreground 
sources, we clearly detect two zsso-dropouts at 3.6pm and 4.5/im, while two others are marginally detected. 
The mid-infrared fluxes strongly support their interpretation as galaxies at z ~ 7, seen when the Universe was 
only 750 Myr old. The IRAC observations allow us for the first time to constrain the rest-frame optical colors, 
stellar masses, and ages of the highest redshift galaxies. Fitting stellar population models to the spectral energy 
distributions, we find photometric redshifts in the range 6.7-7.4, rest-frame colors U — V = 0.2-0.4, V-band 
luminosities L v = 0.6-3 x 10 10 L Q , stellar masses 1 - 10 x 10 9 M Q , stellar ages 50-200 Myr, star formation 
rates up to ~ 25 Mq yr" 1 , and low reddening Ay < 0.4. Overall, the z = 7 galaxies appear substantially less 
massive and evolved than Lyman break galaxies or Distant Red Galaxies at z = 2-3, but fairly similar to 
recently identified systems at z = 5 - 6. The stellar mass density inferred from our z = 7 sample is p* = 1 .6^'| x 
10 6 M Q Mpc" 3 (to 0.3L* =3 ), in apparent agreement with recent cosmological hydrodynamic simulations, but we 
note that incompleteness and sample variance may introduce larger uncertainties. The ages of the two most 
massive galaxies suggest they formed at z > 8, during the era of cosmic reionization, but the star formation rate 
density derived from their stellar masses and ages is not nearly sufficient to reionize the universe. The simplest 
explanation for this deficiency is that lower-mass galaxies beyond our detection limit reionized the universe. 
Subject headings: galaxies: evolution — galaxies: high redshift — infrared: galaxies 



1. INTRODUCTION 

Observations of massive galaxies at high redshift with the 
Hubble Space Telescope and the Spitzer Space Telescope are 
revolutionizing our knowledge of the early formation history 
of stars and galaxies. Blue star forming galaxies at z = 2- 3 
with stellar masses 10 10 -10 u Mq are routinely identified 
from optical imaging (Stei del et al. 1996 a,b,200 4), and have 
been studied in detail (e.g., Papovich. Dickinson. & Ferguson 
1200 ll Shapley et al. 2001,2005), while near-infrared surveys 
at z = 2 - 3 have uncovered substantial numbers of redder, 
more evolved galaxies with larger stellar masses > 1 1 1 M Q 
( FranxetaDEool lYan et alJl2004l iDaddi et alJl2005T) . Some 
of these red galaxies appear to have stellar ages > 1 .5 Gyr, 
im plying that they f ormed most of their stars before z ~ 
5 (Labbeetal. 2005), and suggesting that massive galaxies 
should exist well beyond these redshifts. Direct detection of 
such galaxies would place strong constraints on galaxy for- 
matio n models (e.g.. ISomerville et alJ l2001: Nagamin e et alJ 
2005). 

Tantalizingly, the most recent surveys with the Advanced 
Camera for Surveys (ACS) and the Near-Infrared Camera and 
Multiobject Spectrograph (NICMOS) on the Hubble Space 
Telescope have identified sources at z = 5-6.5 with fairly 
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evolved stars and stellar masses of 1-4 x 10 w M Pl (e.g. , 
lYan et alJ2005llEvles et alJ2005tlpow-Hvgelund et alJ2005l) 
or perhaps more (Mobasher et al. 2005). Critical for these 
results was access to the rest-frame wavelengths longward 
of the B aimer/4000 A break offe red by the InfraRed Array 
Camera (IRAC; Fazio et al. 2004) on Spitzer. Without mid- 
infrared photometry to very faint magnitudes the stellar ages 
and masses of z > 5 galaxies are poorly constrained. 

In this Letter, we extend mass estimates to z = 7 - 8 by an- 
alyzing the mid-infrared fluxes of 6 zsso-dropout candidates 
found in the Hubble Ultra Deep Field (UDF) by Bouwens et 
al. (2004, hereafter B04). These candidates were selected 
from exceptionally deep optical ACS (B eckwith et al. in prep ) 
and near-infrared NICMOS imaging ( Thomp son et alJ 2005). 
and when combined with the ultradeep IRAC data available 
there, offer us an ideal opportunity to verify their reality and 
to study their stellar populations. The stellar masses and ages 
of Z850-dropout galaxies would provide us with the first direct 
look at galaxy formation at z > 7, building on the compre- 
hensive z ~ 6 study (Bouwens et al. 2006b). Where nec- 
essary, we assume an VIm = 0.3, Ha = 0.7, cosmology with 
Hq = 70 km s^Mpc" 1 . Magnitudes are expressed in the AB 
photometric system. 

2. OBSERVATIONS 

The Great Observatories Origins Deep Survey (GOODS; 
Dickinson et al., in prep) observed the UDF with IRAC in 
two epochs, each time integrating for « 23.3 hours in the 
3.6,4.5,5.8, and 8.0/im channels. 6 We estimate limiting 
depths in the combined IRAC images by measuring the effec- 
tive flux variation in random apertures on empty background 

6 This paper uses data release DR3 of epoch 1 
and data relea se DR2 of epoch 2, available from 
|http: / /data .spitzer . caltech . edu/popular /goods /| 
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FIG. 1 . — Images of zgso-dropouts in the HUDFfrom HST/ACS (('775 ,Z85o), NICMOS (/no, #160). and Spitzer/IRAC (3.6/^m,4.5/xm). All sources are undetected 
(< 2a) at ?775 and bluer wavelengths, but we note that zsso-dropout galaxies can have some flux at £775 due to incomplete absorption between Lyman-o at rest- 
frame 1216A and the Lyman limit at 912A. The top two candidates are clearly detected in the IRAC 3.6 and 4.5/xm images, while two others are marginally 
detected. Flux contribution from nearby sources was subtracted in the Spitzer images. Each panel is 4" 1 X 4" 1 in size, or fs21 kpc at z = 7. 



regions. The limits for point sources are 27.7, 27.2, 25.1, 
and 24.9 in the four channels (la, total, 2"5 diameter aper- 
ture). We supplement the observations with deep /f s -band data 
from the Very Large Telescope and Magellan (Labbe et al. 
in prep) and we use an independent reduction of UDF NIC- 
MOS data with improved nois e properties and fewer artifacts 
( Bou wens & Illingworfhl2 "006. hereafter B06). The new NIC- 
MOS images revealed that 2 of the original 6 Z850 -dropouts 
were electronic ghosts of nearby bright stars; hence we re- 
moved them from the sample. 

Matching ACS/NICMOS and IRAC photometry is chal- 
lenging because of the much larger size and extended wings 
of the IRAC point-spread functions (PSFs) resulting in flux 
contamination by nearby foreground sources. Visual inspec- 
tion shows that 2 Z85o-dropouts are substantially blended and 
all 4 are likely to contain at least some flux from nearby ob- 
jects. We have developed a technique to robustly subtract the 
contaminating flux (Labbe et a l., in prep). Briefly, we first 
detect sources with SExtractor (Bertin & Arnouts 1996) in a 
summed NICMOS /iio + #i60 image to determine the light 
distributions at high resolution using the pixels in the "seg- 
mentation" maps. We then convolve these template images 
individually with a carefully constructed kernel to match it 
to the IRAC PSF. Third, we fit all detected sources, includ- 
ing neighbors, simultaneously to the registered IRAC image, 
leaving only the flux scalings as free parameters. Finally, we 



TABLE 1 

Corrected Photometry of z 850 -dropouts in the UDF. 



ID 



~S5(! 



#160 



3.6/xm 4.5/^m 



964 > 29.6 26.8±0.2 26.8±0.2 27.0±0.9 26.1±0.2 26.4±0.5 

1417 28.0±0.2 26.7±0.2 26.3±0.2 26.1±0.3 25.3±0.1 25.5±0.2 

950 > 29.6 26.9±0.2 26.6±0.2 > 27.2 27.0±0.5 26.4±0.5 

1125 29.2±0.7 27.3±0.4 27.2±0.3 > 27.2 26.9±0.4 > 27.2 

NOTE. — Object IDs are from B04. All magnitudes are in the AB system. 
Optical/near-infrared fluxes were measured in 0"9 diameter apertures and 
IRAC fluxes were measured in 2"5 diameter apertures. We corrected the 
fluxes for missing light outside the aperture assuming stellar profiles. The 
corrections are 5 — 20% in the optical/near-infrared and a factor 1.75, 1.81 
for IRAC 3.6/im, 4.5/jm, respectively. We a djusted the NICMOS Ji in and 
#160 zeropoints by -0.16 and -0.04 following de Jons et al. 12006). In the 
improved NICMOS reduction (B06) object 1417 falls within the Z850-dropout 
selection criteria, whereas in B04 it fell just outside. The upper limits are la. 
None of the galaxies are detected in V^, (775, and 5.8/-im to la limits of 30.6, 
30.5, and 25.1, respectively. 



subtract the best-fit images of all neighboring sources. 

After cleaning the IRAC images, we performed con- 
ventional aperture photometry in 3.6,4.5,5.8, and 8.0/im 
bands in 2"5 diameter apertures. Photometry in the ACS 
fi435,V606, '775,2850, NICMOS /no, #160, and K s bands was 
done in 0."9 diameter apertures and we obtained magnitudes 
and limits consistent with B04. We summarize the photome- 
try in Table 1 . 
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FIG. 2. — Examples of the deblending procedure in the IRAC 3.6/xm im- 
age as performed separately for the first (top row) and second (bottom row) 
epoch images from GOODS. We fit the original images (left) using the high 
resolution NICMOS data to construct models for the nearby neighbors (mid- 
dle) which are subtracted resulting in the cleaned images (right). We re- 
cover consistent fluxes for both epochs after subtraction. The panel sizes are 
4"lx4"l. 



3. MID-INFRARED FLUXES OF Z 850 -DROPOUT SOURCES 

Figure 1 shows the i/S77ACS+NICMOS and the combined 
Spitzer/IRAC images of the z&so dropouts. Two objects (ID 
964 and 1417) are unambiguously detected in 3.6/im (5-8<r) 
and in the slightly shallower 4.5/im (2-3er). Two others (ID 
950 and 1 125) are only marginally detected, but are probably 
real as the sum of their 3.6 and 4.5/im images reveals a visi- 
ble source. Unfortunately, the IRAC observations are not deep 
enough to definitively confirm or reject the reality of the un- 
detected sources. None of the candidates are visible at 5.8 and 
8.0/im. To evaluate the robustness of the deblending photom- 
etry, we performed the procedure independently on the first 
and second epoch IRAC data (see Figure 2). Reassuringly, 
we measure consistent fluxes and we detect the brightest, most 
promising sources 964 and 1417 in each dataset individually. 

The 3.6/im magnitudes are faint, ranging from 25.3 to 27.0, 
with Hi(,o- 3.6/im colors in the range -0.4 to 1.0. The IRAC 
detected objects 964 and 1417 are the reddest, showing a fac- 
tor of ~ 2 rise in f v flux densities between H\qo and 3.6/im, 
while the SEDs are flatter at 1.1— 2.2/zm and 3.6 — 4.5/im 
(see Table 1). The rise at 3.6/im is similar in strength to 
what is found in spectroscopically confirmed galaxies at z = 6 
jYan et al.l2005llEvles et al.l200 5). and suggests the presence 
of a substantial redshifted Balmer break, indicative of evolved 
stellar populations. 

4. STELLAR POPULATIONS OF Z « 7 GALAXIES 

We fit stella r pop ulation synthesis models of 

Bruzual & Chariot (2003) and a iCalzelti et.alJ (2000) 
obscuration law to the broadband fluxes to constrain the 
stellar populations. The models assume Solar metallicity 
and a Salpeter initial mass function (IMF) between 0.1 and 
100 Mq. We explore three different star formation histories 
(SFHs): a single age burst (SSP), an exponentially declining 
star formation rate (SFR) with a timescale of 100 Myr (tioo), 
and constant star formation (CSF). We leave the redshift, 
mass, age, and extinction as free parameters. 

We find acceptable fits for all sources (see Table 2 and Fig- 
ure 3) and obtain confidence intervals on the parameters with 
Monte Carlo simulations. The best-fit redshifts vary from 6.7 
to 7.4 and most Monte Carlo solutions are in a narrow range 



TABLE 2 
Best-fit Model Parameters 



Solar metallicity 



ID 




Mass 


Age,,. 


Av 


SFR 








(1O 9 M ) (Myr) 




(Mq yr 1 ) 








Single 


Age Burst (SSP) 






964 
1417 


7 9 +0.1 
'• z -0.2 
£ 7 +0-2 
6 - 7 -0.2 


1-8 


40 
81 


0.0 
0.1 


0.0 
0.0 


1.5 
0.8 


950 


7 7 +0.1 
'■ z -0.1 


1 o +0 - 2 

l.U _q i 


25 


0.0 


0.0 


1.4 


1125 


7 +0 2 

' -0.3 


7 +0 ' 4 
u -' -0.1 


32 


0.0 


0.0 


0.9 


Exponentially Declining SFR (e-folding time r 


= 100 Myr; tioo) 




964 


7 7 +0.1 
'■ z -0.1 


3.4 *>;« 


133 


0.0 


5.2 


1.2 


1417 


£ Q +0.1 

68 -0.2 


9.3 «j 


176 


0.2 


10.7 


0.7 


950 


7 T +0.1 

'• 3 -0.2 


0.6 _ 3 


35 


0.0 


9.4 


1.1 


1125 


7 1 +al 

' • 1 -0.2 


9 +o s 


77 


0.0 


4.4 


0.8 


Constant Star Formation (CSF) 


964 


7.4 1i 


4.2i|f 


320 


0.0 


7.3 


1.3 


1417 


6.8^ 


11.5 31 


254 


0.4 


25.5 


0.9 


950 


7 o +0.1 
'• 3 -0.2 


7 +0 <i 
U. / _ 4 


40 


0,0 


9.9 


1.1 


1125 


7 - -O2 


j q+1.3 
lM -0.6 


127 


0.0 


4.5 


0.8 



NOTE. — The best-fit paramete rs for Bruzual & Chariot (2003) stellar population 
models and a Calzetti et al. 2000) obscuration law. We omit the B and 8/jm band 
from the fit as they are not deep enough to constrain the models, leaving 9 filters 
(^606, '775, 2850, 7no,77i60,^,, 3.6/im, 4.5/im, and 5.8/^m), and 4 free parameters (red- 
shift, stellar mass, age, and visual extinction Ay). We assume Solar metallicity and a 
Salpeter IMF from 0. 1 - 100 Mq ; adopting a Chabrier IMF would result in similar ages 
but ~ 1.7x lower stellar masses. Sub-Solar (1/50 Zq) metallicities result in 30% higher 
ages and 20% higher masses. We consider redshifts from to 10 and Ay from to 2. We 
explore 3 star formation histories (SFHs): SSP, Tioo, and CSF. The 68% confidence in- 
tervals are obtained from Monte-Carlo simulations and the reduced x^j values assume 5 
degrees of freedom. We impose a minimum best-fit age of 25Myr to avoid unrealistically 
young solutions. Finally, we weight the best-fit age with the SFH to better represent the 
age of the stars comprising the bulk of the stellar mass (Age,,) For CSF the correction is 
0.5, for exponentially declining S FHs it is [age — r + Texp(-age /r)]/[\ -exp(—age/r)], 
while SSP requires no correction Forster Schreiber et al. 2004). 



around the best fit indicating that the redshift is well con- 
strained. Old stellar populations at z ~ 1 fit the data poorly 
as they do not reproduce the strong break across the 1775,2850, 
and Ji 10 bands and the blue near-infrared continuum. Even so, 
we note that 3% and 1 1% of the solutions for object 1417 and 
1 125, allowed a redshift of z ~ 1 when the random flux vari- 
ations "conspired" to suppress the break. The best-fit models 
have average rest-frame optical colors of U—V = 0.4 (0.2) and 
V-band luminosities L v = 2.3(1.0) x 10 10 L o for the IRAC- 
detected (undetected) sources. 

Determinations of the stellar population age and mass de- 
pend on the assumed star formation history and metallicity. 
For the whole sample, tioo models fit the best, with ages 
of 40-180 Myr and instantaneous SFRs of 4 - 11 M yr" 1 . 
Converting the rest-frame 1500 A luminosities directly into 
(absorption-corrected) SFRs (Ma dau et al.lll998T) results in 
similar values. Nevertheless, the degeneracy between age and 
dust prevents us from placing firm limits on the ages and 
SFRs. The stellar masses are generally better constrained. 
As expected, the IRAC-detected galaxies 964 (z p h = 7.4) and 
1417 (z p h = 6.8) are the reddest, most massive, and oldest in 
the sample. The average uncertainties on the masses are ap- 
proximately a factor of 2 - 3. 

The extreme SSP and CSF models set lower and upper 
boundaries to the masses, ages, and SFRs, whereas assuming 
sub-Solar metallicities (1/50 Zq) instead of Solar results in 
30% higher ages and 20% higher masses. Because all SFHs 
and metallicities provide equally acceptable fits to the data, 
we will hereafter adopt the mean of the SSP and CSF models 
and both metallicities as our fiducial values. We then find ages 
of 50-200 Myr, masses of 1 - 10 x 10 9 M Q , low reddening 
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FIG. 3. — The observed SEDs of the IRAC-detected z ~7 galaxies and 
the best-fit Solar metallicity models for each star formation history. The ver- 
tical and horizontal bars mark the Icr flux uncertainty and the width of the 
passband. The arrows are lcr upper limits. The sources show evidence for a 
break between the near-infrared and mid-infrared fluxes, which is well fit by 
a fairly evolved stellar population. The fits generally require little reddening, 
although galaxy 1417 also allows for moderate amounts of reddening and 
fairly large star formation rates (up to Ay = 0.4 and SFR ~ 25 Mq yr _1 ), 



Ay < 0.4, and star formation rates of 3 - 12 M Q yr . 

5. DISCUSSION 

Using the GOODS dataset (Dickinson et al. in prep), we 
have estimated and analyzed the Spitzer/IRAC mid-infrared 
fluxes of 4 zggo-dropouts candidates, which were identified in 
the UDF by B04 and remeasured more accurately in B06. The 
sources are rare, with a surface density of 0.7 arcmin" 2 , and 
very faint, with observed magnitudes of Hj6o = 26-27 and 
3.6/i = 25-27, placing them well beyond the spectroscopic 
capabilities of current telescopes, but in reach of future facil- 
ities such as JWST and ALMA. IRAC directly confirms the 
reality of two sources, while two others are marginally de- 
tected. Modeling of the broadband fluxes strongly supports 
their interpretation as z ~ 7 galaxies with substantial stellar 
masses 1 - 10 x 10 9 and ages 50-200 Myr. 

Using the redshift selection function for the z^so-dropout 
sample (B04), we obtain an effective volume of 9000 Mpc 3 to 
0.3L* =3 and we infer a stellar mass density of = x 



10 6 M Q Mpc~ 3 . Comparing to the stellar mass density at 
lower redshifts, computed using similar techniques to sim- 
ilar luminosities, we find a continuing decrease where the 
densi t y at z ~ 7 is 95% , 22% of that at z = 6,5 JYan et alJ 
120061 Istark et alJl2006l) . Recent Smooth Particle Hydrody- 
namics (SPH) simulations in a ACDM universe predict stel- 
lar mass- and number densities for massive (> 1.8 x 10 9 M Q ) 
galaxies of 0.9 x 10 6 Mp, Mpc " 3 and 2.5 x 10" 4 Mpc" 3 (SPH 
G6 run, Naaamine et al. 2005), remarkably close to the min- 
imum estimates for our sample (0.7 x 10 6 M Q Mpc" 3 and 
2.2 x 10~ 4 Mpc" 3 to the same mass limit). However, we 
cannot exclude additional stellar mass residing in massive 
non-starforming or dust-enshrouded galaxies, which the Z85o- 
dropout criteria would have missed. 

The ages of the two most massive (IRAC-detected) z%sq- 
dropouts suggest that the bulk of their stellar mass formed 
at even higher reds hifts z > 8, during the epoch of cosmic 
reionization (Spergel 'et all2006l) . We can place a simple up- 
per limit on the contribution of high-mass galaxies to reion- 
ization, by calculating the maximum SFR densities implied 
by the observed stellar masses and ages for these objects. 
Taking the maximum masses and dividing it by the mini- 
mum ages (~ 60Myr), we infer a substantial SFR density 
0.04 M Q yr _1 Mpc" 3 . This is higher than calculated directly 
from t he rest-frame UV lumi nosities at z = 6 - 7 to the same 
limits (Bouwen s et alJ200 6b). but still more than 3 times too 
small to rei onize the universe a t the lowest probable redshift 
Zreion = 8.6 (Spergeletal . 2006) for cano nical assumptions 
( Mada uet all 1 9 99 : Bouwens et al. 2006b, see also Yan et al. 
2006). 

The simplest explanation for the lack of ionizing pho- 
tons provided by high-mass galaxies at z = 7 is that low- 
mass galaxies beyond our current detection limits were pri- 
marily responsible for reionization t Lehnert & Bremerll2003l 
Yan & Windhorst! 120041 iBouwens et alJl2006bl) . Other pos- 
sibilities are that massive galaxies are missing from current 
surveys and are unaccounted for in the models, or that the ob- 
served galaxies had a top heavy IMF, which would increase 
the ionizing efficiency per unit stellar mass. Obviously, in- 
completeness, sample variance, and large scale structure dom- 
inate the uncertainties in our results. Larger, very deep near- 
infrared surveys from the ground and from space will address 
this issue in more detail. 

We are most grateful for the efforts of the GOODS team 
who provided the data used in this analysis. We thank Ken 
Nagamine for providing his numerical simulations and the 
referee for helpful comments. IL is supported by a fellowship 
from the Carnegie Institution of Washington. RJB and GDI 
acknowledge support from NASA grant HST-GO09803.05-A 
and NAG5-7697. 
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